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The interactions of gravitational waves with interstellar matter, dealing with resonant wave- 
particle and wave-wave interactions, are considered on the basis of magnetic-type Maxwell- Vlasov 
equations. It is found that the behavior of the fields, involving the "gravitoelectromagnetic" or " 
GEM " fields, the perturbed density field and self-generated gravitomagnetic field with low frequency, 
can be described by the nonlinear coupling equations Eqs. 16.10H - 16.12ll . Numerical results show 
that they may collapse. In other words, due to self-condensing, a stronger GME fields could be 
produced; and they could appear as the gravitational waves with high energy reaching on Earth. In 
this case, Weber results, perhaps, are acceptable. 
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I. INTRODUCTION 

It has often been noted that magnetism can be understood as the consequence of electrostatics plus Lorentz 
invariance. Similarly, Newtonian gravity together with Lorentz invariance in a consistent way must include agravit- 
omagnetic field. This is the case of gravitoelectromagnetic form of the Einstein equations in a medium P, 0, HE! 0- 
It is well known that in the slow-motion limit of general Orelativity, accurate to post-Newtonian (PN) order, where 
the field equations can be reduced as [||: 

V 2 </> = AnGp, (1.1) 

V 2 A = 167rGpv/c; (1.2) 
moreover the harmonic gauge condition and the force on a unit mass reduced to 

4 d 

V-A+ = 0, 
c at 

dv _ , 1 dA v _ 

~_V0 — + - x Vx A, 

at c ot c 

where <f> is the Newtonian gravitational potential, Aj/c = gio is the mixed metric and G is in the Newton's constant. 
One defines the GEM fields via 



1 8 A 

E 9 = -V^-~ — , VxA = B 9 (1.3) 



in direct analogy with electromagnetism; it follows from these definitions the Maxwell-type field equations in the 
continuous medium0,Q 

VxBs = T js ^r V-D 9 =47r/> s , (1.4) 
1 r)R 

VxE 9 = V-B 9 =0 (1.5) 

c at 
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where p g = —Gp is "matter density" and j g = — Gpv = p g v is "matter current". Then force on a unit mass has 
Lorentz-type form 



dv _ v 

at c 



E 9 + : xB 9 . (1.6) 



And all terms of 0(u 4 /c 4 ) are neglected in the above analysis. Now one decompose the total matter density and 
current in Eq. (|1.4|) : p g = p g + p g o (or p = p + po), j g = j g + j g o, in which the p g o(ov po) and j ff o are the parts of 
a external field source in the medium, correspond to a local additional mass disturbance, for example, by nonlinear 
interactions between the fields and the medium, or some test particles. In this case it is convenient to definite the 
vector quantity D s (£, r) through ( in what follows one neglects the marks " ~" in p g and j 9 for simplicity ) 



D 9 = E s + 4tt 



J dt%(t',r); (1.7) 

— oo 

by use of Eq. I|1.7fl and the mass continuity equation, 



— + V-p Q v = 0, (a =5, g ) 



Eas. H1.4|) and (|1.5|l are deduced to 

16tt . 4 9D 



V x B g = — j g0 + -^, V-D s =47rp g0 , (1.8) 



1 r)"R 

VxE 3 = V-B s =0. (1.9) 



A dimensional analysis of equations (|l.lf) . (|1.3I) . I|1.7|l and l|1.9|l suggests that it is insightful to use the CGS units, 
since these field variables have the following dimensions: 



thus, one has 



by substituting 



[D g ] ~ [E g ] ~ [B ff ] ~ [V0] ~ v / a , 



[l£/G] ~ [B g /G] ~ ~ M ; 



E = -E S /4VG, B = -B g /WG, D = -D 3 /4VG, (1.10) 
j = -j 3 /Vg = VGpv, p = -p g /VG = VGp, 

the GEM field equations can be expressed in the standard form, 



4tt ldD ldB , 

VXB -T J0+ C^' VXE --C^' (L11) 



V-D = 7r/5 , VB = 0; 



(1.12) 
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then Eas. (|1.7|) and l|1.6() become 



D = E + 7T / dt'j(t',r) 



(1.13) 



dt 



E + - x B 

c 



(1.14) 



To consider the responses of the medium on the GEM fields, one must introduce material relation, which describes 
the GEM properties of the medium. In view of Ea. tjl.lHfl . the states of the medium not only depend on a given time- 
space point(£, r), but also depend on previous times and at any point of the medium. Hence, by general reasoning in 
physics ( independent of a specific model for the medium ) one can state that this is a no-local linear relation in the 
limit of linear response ; whose Fourier representation is 



ji(ui,k) = aij(u>,k)Ej(w,k). (1.15) 

In special, for common continuous medium, where the "spatial dispersion" (dependency on k) is not important , the 
relation l|1.15fl is deduced to (<Xy — > <7<5y ) 

= cr(w)E(w), 

i.e. an "Ohm's gravitational law" , where a is gravitational conductivity (see Refs [10,11] for details) . Then by 
use of Eqs. I|1.13|l and (|1.15|) . one has 

U i (cj,k)=e ij {u,k)E j (ij,k), (1.16) 
here ey, called the dielectric tensor, is determined as 



e«(w,k) =<% + — <7«(w,k) {cj^O). (1.17) 

In the next section we shall show that the material relations Eas. (|1.15fl and (|1.16|) are relevant for GEM fields. 

It has often been noted that the effects of GR can largely be understood and treated easily within the GEM 
framework. For example, the famous Lense-Thirring precession effect of GR is simpler and clearer to use the GEM 
equations [13 ; and the gravitomagnetism is a useful insight for understanding the Schiff effect |l3J and the "Faraday" 
effect fl4| . In special, just from our electromagnetic experience, we can infer that at a distance should present grav- 
itational waves predicted by GR. Therefore, the existence of the GEM fields is equivalent of the existence of the 
gravitational waves predicted by GR|13J. 

The many efforts that have been made to detect gravitational waves have so far given no convincing evidence 
that they have actually been seen. In the late 1960s and early 1970s, Weber announced that he had recorded 
simultaneous oscillations in detectors 1000 km apart, waves he believed originated from an astrophysical event. But 
many physicists were suspicious of the results that were several orders of magnitude higher than were theoretically 
predicted. In addition, observed high energies of gravitational waves by Weber have not been confirmed by these 
independent measurements. This is perhaps, due to the fact that gravitational waves with high energies are very rarer 
than many physicists had expected. On the other hand, it is highly possible that the gravitational waves can get 
largely increase rate as the waves interact with interstellar matter, so that Weber's result, perhaps, are acceptable. 

In order to study the nonlinear effects on a very larger scale, where the mean free path for collision between the 
particles of the matter is larger compared with the characteristic length appearing in the problem, a kinetic treatment 
is required: The physical systems should be described by the magnetic-type Maxwell- Vlasov equations. The small 
amplitude and high-frequency approximation are employed for the gravitational waves, such that perturbed techniques 
can be applied. The effects deal with resonant wave-particle and wave-wave interactions, inhomogeneties of the matter 
distribution and nonlinear self-collapsing. A previous work|15j is just devoted to the study of the nonlinear interactions 
in the absence of gravitomagnetic field. The present paper considering gravitomagnetic component is a generalization 
of the work. And our another work^^ is its analog of electromagnetism. 
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In Section II we establish the description of collisionless kinetics and give the linear effects. Then we present the 
nonlinear equations of the fields with low frequency and high frequency, starting from GME equations in Sections III 
and IV . In Section V we focus on the motion of matter disturbed by GME fields. As a result, the nonlinear controlling 
equations, involving the coupling of GEM fields with the perturbed density field and self-generated gravitomagnetic 
field, are presented in Section VI . A numerical integral of the controlling equations is given in Section VII. Finally 
we sketch the conclusions, stressing the possibility of detecting gravitational waves. 

II. KINETIC DESCRIPTION AND LINEAR EFFECTS 

Taking into account the presence of massive dark matter in the universe, one should treat a two-component self- 
gravitating system . The collisionless Boltzman equations for the distribution function f^are[l5| 

^+v^ + (a + F).^=0 (. = 1,2), (2.1) 
where a is the non-gravitational term, F is Lorentz-type force: 



E+ - x B 

c 



(2.2) 



In which gravitoelectric field E and gravitomagnetic field B satisfy Eas. (jl.llfl and H1.12|l . and the "charge" is 

q a = -4:VGm a . 

The density n a and the current density j are connected with the particle distribution through 

n»(r,t) = J U*,v,t)j0-, (2.3) 

and 

jM) = -|£ / fcv/«(r,v,<)^3. (2.4) 
We can assume to the gravitation effects 

p 2 = n Q m 2 < p 2 = n Q m 2 , (2.5) 

where p\ denotes the density of bright matter, and p 2 the density of dark matter. We divide f a , E and B into two 
parts: unperturbed and perturbed parts, 

fa = f* + fl, E = E fl + E T , B = B i? + B T . (2.6) 
As the f a is closely coupling with GEM, the perturbed distribution can be expanded in power 

/J = £/ TW (2-7) 

i 

provided that the perturbed field E T is weak, i.e. 

IE T I 2 

W = { 7 F «h (2.8) 
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where the index i indicates the i-th power of E T . In this case that equation for the unperturbed state 

f)fR a f R a fR 

"Jot , v .^ + (a + F fl).^£L =0 (2.9) 



reduced to 



dt dr dp 



QfR 

= 
dt ' 



with a relevant solution 



tR dp y dp z _ (2tt 



,1/2 



ti*. - / /<?7^ = T^^^oe"^. (2.10) 
:P 1 \2-kY (m a v Ta ) 



Substituting Ea. (|2.6|l into Ea. (|2.1jl and subtracting Ea. <|2.9jl . putting a + g R ; = 0, yield 



3/1 , v ,^ +f t ,dfl +F T ,dfl^ Q (2 n) 



Si <9r dp dp 

Substituting Ea. (|2.7ll into Ea. (|2.11|l and expanding A = (F T , / Q )in a Fourier series 



A(r, v, t) = J A k e- tuJt+ *- r dk, (2.12) 



we get from Ea. H2.11|) : 



Ak = Ak u , dk — dkdw, 



iu< kY)l:[',' ) =lFl-^^S(k-k 1 -k 2 )dk 1 dk2, (2.13) 



g f T(l) 

i(u - k- v)/^ = / F£ • -^5(k -h- k 2 )dhdk 2 , (2.14) 



Dividing j into to two parts, 



i(w - k • v)/Jjf } = y • -^^5(k -ki- k 2 )dk x dk 2 . (2.15) 



j=j fl +j T , (2-16) 



and expanding j T in powers of -E T 

jno__» E | tTj? «* (, 17) 



(2t) 

one can obtain the linear current from Eas. i|2.13|) . H2.10|) . and l|2.17|l 



3k,i = °y( w > k )^j ' 



(2.18) 
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This is the relation Eg . (|1 . 15|) : where 



flW 1 V rViq " 



0/* dp 



i(co — k • v + ie) dp s (2tt) 3 



(2.19) 



and here we have into account the slow change in f^, i.e. f^ k2 ~ fa^i^z) the term ie arises from the Landau rule. 
Similarly, we get the nonlinear currents (up to the third order) from Eas. l|2.14l) and (|2.15|l 



; J Sl k2 k3 ElEldk lC lk 2 6(k -kt- k 2 ), 



(2.20) 



*k,ki ,k2 



Hi E k E kJ ( k - fc i -ks-k 3 ) dkidk 2 dk 3 



(2.21) 



with 



1 



(®ki ' 9 9 p) ( e k 2 ' 9p) K dp 



Sfe, fcl , fc2 / {u J k . v + ^) 1 k2 . v | fe ) /a (27r) 3 > 



(2.22) 



Get 



1 



v • dp 



1 



a 



i 



1 '•'./ (2tt) 3 (w-k-v + ie) V k ' dp J [(u-wi)- (k-ki) • v + ie] 
9 \ 1 / CT d \ R 



; k 2 ' 



dp J (lu3 — k 3 • v + ie) \ k3 dp j 
On the other hand, we can obtain the field equation from the Maxwell-type GEM equations 



/ n>2 



where 



(2.23) 



(2.24) 



4 = C k = k)e^er ? + — (k • e£)(k • e£* 



(2.25) 



is "dielectric constant" for <r mode (longitudinal mode, a = I and transverse a = t) and E fe " = S fc "e£ with e£ -e?*, = 1. 
Using Eas. ((TT7jl . (EHJand l|ZTU|l . one has 



4 = e ^.(w,k) 



fc 2 



1 ^ fc 2 



j 2 

pa 



"i'Q 



V \/2kvTa 



(2.26) 



1 wf 



^ A; 2 LO \y/2kv T a 



LO 



(2.27) 



with 



(2.28) 
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where 



is dispersion function pl|: 



\V2kvTi 



X/y/TT 

x — £ + is 



^2 



(2.29) 



13 2 

Z(x) Rill — r + — r - iypnxe~ x , X > 1, 
2x 2 4a; 4 



(2.30) 



For high-frequency field 



Z(x) w 2x 2 — iy/irxe x ~ 2x 2 — iy/nx, 



(2.31) 



i.e.x 3> 1, therefore, neglecting the damping term, one obtain from Eas. (l2.2ti|l and l|2.27|l 



(2.32) 



4 = i 1 



W„9 w. 



2 



0^ or 



and 



(2.33) 



For the low-frequency fields, the following conditions are met 



2 ' 



(2.34) 



v T2 > ui'/k' <C vti , w'-^Wpi, 



Then 



(2.35) 



4' 



i + - 1) + i) 



U T1 



' J p2 



U T2 



(2.36) 



<>2 



^ 1 7T W p 2 Wp2 



2 u/ fc'VT2 



(2.37) 



Taking Ea. l|2.34|l into consideration, one gets the linear dispersion relationship from Ea. H2.24|) for the transverse 
oscillation with high-frequency: 



U, 2 = k 2 C 2 - (i0 2 pl + LU 2 p2 ) 



(2.38) 



Usually, there is not steady longitudinal high-frequency mode on basis of Eq. (|2.33l) , which consistent with the analyses 
for the electro-gravitation kinetics [l^. Consider that o nly the oscillations with frequency is close to the proper 
frequency of the medium, which is similar to plasma case jl7l| : this means 



2 2 i J"2 2 2 
U) ~UJ pn + k C ~UJ p2 , 



2 2 i 2 T 2 2 

Upn = ^pl + ^p2 > k C 



(2.39) 



Ea. (|2.39|l is a branch of dispersion relationship [Ea. (|2.38|l ]. In what follows we omit the mark '~' in Ea. (|2.39|l for 
simplicity. 
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III. LOW-FREQUENCY TRANSVERSE FIELD EQUATION 

First of all, let us study the case of the low- frequency transverse field — E^* = E^ s in Ea. (|2.24|l . For the 
second order nonlinear current, because of delta function in Ea. (|2.2U[l . it must met k = ki + k 2 , u = oj\ + ui 2 ; since 
(k, (J) belong to low- frequency wave, then lo\,uj 2 must be high-frequency and are of opposite sign. In this case, the 



second order term in Ea. H2.2(J[) becomes: 



E T(+) E T{-) E T(-) E T(+) 



where upper indices "+" and "-" denote the 



positive and negative frequency parts of high-frequency perturbations, respectively. Hence 

ti ] = E / s tSlM + Hk - h - k 2 ) dk ld k 2 , (3.1) 



where 



s «(f) _ 1 3 f , r e ki ' dp e L ' dp fa dp 

3 fe,fc 1 ,fc 2 - ,9a / V 3 . {6.Z) 



4 a J (w - k • v + is) uj 2 - k 2 • v + is (2tt 
Using the substitution k\ — > k 2 , £; 2 — * fci in the second integral term of Eq. I|3.1|) , it yields 

i 2) = E / ( S "L + S?2, fcl )< C+ X ( " , *(* ~ fc i - k2)dk ld k 2 ; (3.3) 

a J 

substituting Ea.lpT3)) into Ea. $Z3ty yields 

(fc 2 c 2 - W 2 4)^ TS =4^^ | ^^^-^(fc-Ai-feJdfcjdfca, (3.4) 



where 



i ^ e k 2 • gjj | _ t 9 j'l 1 / tt % 
kl 3p w 2 - k 2 ■ v + is k2 dp u>i - ki • v + is J (2tt) 3 ' 



with 



v ,, k v ., . v • C> L 



eL = (l-— -)«£ + (— ^)k; ,:;.(,) 
and the second order distribution function is 

fal = J ^A,*X (+)i C ( ~ )<J ( fc - fc i - k 2 )dk 1( ik 2 , (3.7) 

where 

y a _ 2 1 f-t 9 e L 2 ' ^ f 9 1 f fl /ooN 

^-k-v + ze 6kl '^a; 2 -k 2 -v + ie + ^"^wi-kx-v + ie /q ' ( ' 



9 



IV. HIGH-FREQUENCY TRANSVERSE FIELD EQUATION 

If Ef? in the left-hand side of Ea. (|2.24fl is high-frequency field, = E^ + ', the quadratic terms in the Ea. H2.20|) 
should be the product of high-frequency and low- frequency fields:. E% E k — E k h E k s + E k s E k h . The three- 

.(3) 

field product included in the current y k can be expressed in terms of high-frequency fields and cubic a mixed 
product of high-frequency and low- frequency fields. Using Ea. l|3.4|l . then, in fact, this mixed term is the prod- 
uct of four high-frequency fields, which is higher order . Then E k E^E^ « E k E k E k . Due to the factor 

[(a; — uji) — (k — ki) ■ v + ie]" 1 in the Ea. l|2.23[) . its contribution to the jjj, 3 ' is important if E k h is the positive high- 
frequency fields. As a result 



p T p T p T _ p T(+) 



E T k ^El^ + El^E T k ^ 

fc 2 k 3 ' k 2 k 3 



Therefore, we gets similarly the high-frequency field equation as follow 



(fc 2 c 2 - w 2 4) El {+) = iniuj 



e fsitkA^K 36 ( fc - fc i - ^ dk ^+ 



E / Gt% MM E T k ^E T h ^E T k ^8 (k - fcr - fe 2 - fc 3 ) dMfcd** 



(4.1) 



where 



fe,fcl,fc2,fc3 




and the expression forS 1 ^^ fca is the same as fe 2 (see Eq. (|3.5^ . except that is low- frequency now. And E k s is 

low- frequency fields, but it may different from the low- frequency fields E k s in the left-hand of Ea. (|3.4|) . According 
to semiclassical theory, the fusion and decay interactions in Eas. i|3.4|) and (|4.1|l can determine the field intensity with 
low- frequency, N k s ~ \E k s \ 2 , in other words, they differ by a phase factor e 1 ® . 

The symbol ^ in those expressions above implies adding the contribution of dark matter(a = 2) and bright 

matter(a = 1). From Eas. l|3.5l and 14.21) one can see that /^rfp ~ n , q a ~ m Q , S a oc m a , G a oc m a , namely that the 
matrix elements of interaction are proportion to the mass of particles. As the assumption of the sark matter mass 
far larger than the bright one[see Ea. H2.5P ]. we can neglect the contributions of the bright matter. Therefore we write 

c«W o«M fln H r a{t) ^ 2(t) ? 2(t) »mc\ r 2(t) 

°k,k 1 ,k 2 ' k,k 1 ,k 2 < ~ r k,k lt k 2 ,k 3 °k,k 1 ,k 2 ' J k,k 1 ,k 2 duu Lj fe,fci,fe 2 ,fe3 ' 

In order to get the field equation in spectrum space for nonlinear interaction up to the third order, we must estimate 
in detail the integral value of matrix elements^ ^ fc2 and Grj^ k2 k3 in Ea. l|4.1|l . Integrating Ea. l|4.2|l by parts and by 
using of Ea. (|2.39[) and 1)3. 8p. we get 



^G 2 k % k2 k3 « 4^ 2 2 e- • J ^-Imm^ <4.3) 



and for k2 , in which u> and Wi are high-frequency and L02 is low-frequency, after integrating by parts, it reduces 



to 



4 ww 1 w 2 m 

Then first integral of Eq. I|4.1|) becomes 



- 177^^ • K >< >< <)] ■ (4-4) 
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7772 C 

where 



Airiuj J Sl^ iM E T k ^ +) El 2 s 5 {k - h - k 2 ) dk x dk 2 (4.5) 

5 {k — k\ — k 2 ) dkidk 2 . 



% a X B k 2 



B t 2 = — X K S (4-6) 

L0 2 

is low-frequency magnetic filed, which are produced by the fields with positive and negative high-frequencies. And 
using Ea. H3.7f) . the second integral of Eq. (|4.1|l becomes 

(2) 

,2 / „** ?1?T(+) n k-k 1 



4^ 2 /e--(E^^)dfc l7 (4.7) 

77-0 



(2) . 

where n k , is the second order of perturbed density, 



Therefore Ea. (|4.1|) is reduced to 

,(2) 



(JbV - cohi) E^ + ) = ^ p2 [ (^^) dkl + 4^ p2 / (e£< + > x B^)*. (4-9) 

J ^0 777-2 C J ^ ' 

As a result, we obtain from Ea. 14.9f) 

2i <9E(r,i) c 2 . N , 

v ' ; —V x V x E(r,t) (4.10) 



LO p2 dt uj p2 
= -2E(r, t) + 4 n(2)(r '^ E(r, t) + 4 tg2 E(r, f) x B s (r, i) , 

77 m 2 CLU p2 

here E(r, t) is the envelope for the high-frequency fields, 

E(r,i)e^ 2t = I *E^ (+) e-' iwt+4k - r dfc, (4.11) 



and because of slow change in E(r,i), we have neglected the term (J^-)E(r, t). 



V. PERTURBED DENSITIES 



It is possible that two transverse fields with high frequencies can produce a longitudinal low-frequency field. In this 
case, the density perturbation of the first order can be exited. Then according to Eq. (|3.4f) . one has 

a J 

~ 2 f / 1 — 2f £ 1 

where the coupling matrix S k k k2 is the same as S k k k2 when e£* — > k/fc. the density perturbation produced by the 
longitudinal low-frequency wave is 

,(!) 



T(i) dp _ q-2 E TS(l) f 1 k ' IF' rf P 
(2tt) 3 ~ * fc J ku-k-v + ie (2tt) 
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Taking account of the integral expression of ej^fsee Ea. (|2.26(l ]. we obtain 



g2 k 
i 7tq| 



E 



TS(l) 



(5.3) 



Using (e^ — 1 ) = e l k — and e l k E k b{l> , as the first-order approximation, is zero[see Ea. (|5.1|l ]. then 



J .1(0 



.1 P TS(1) 



= g2 fc c l(l) E TS(l) 

k i irq% k k 



(5.4) 



Similarly, one estimate the matrix as follow: 



52(0 

J k,k 1 ,k 2 



kuJ ( 2(0 Aft t ^ 



(5.5) 



So the total perturbed density produced by two transverse fields with positive and negative high-frequencies, is 



(1) , (2) 
n k + n k 



(5.6) 



no _l j 



= -4- 



7rm2o;, 



p2 



- - l) / • E I 2 ( ~^ (* - fci - fca) dfcidfe. 



Then, in the coordinate representation, Eq. l|5.6|l becomes 



n'(r,t) = 4 



7rT 2 



(5.7) 



In the case of taking account of longitudinal low-frequency fields, the first term to the right-hand side of Eq. 14.9(1 
will be added to the coupling term of longitudinal fields, 



4™ / Sl% M El i+) El S{l) 6 (k-h- k 2 ) dhdk 2 



(5.8) 



where the coupling matrix S^2 1 k 2 ' is s i muar to fc 2 ( e k 2 — * if ) m E( l- *4.5|) . its estimated value is 



?2(0 



1 , ret-e* 

-A 



. JL f 

k 2 Op*' 2 



dp 

wni2 ^2 - k 2 • v + ie (2tt) 3 



And Eq. (|5.8|l is reduced to 



(5.9) 



T{+) 



,T(1) 

(2tt) 



E 



„(!) 

r(+) "fc-fc 



fei 



"o 



(5.10) 



Taking this term into account, one gets from Ea. l|4.9|l 



2t SE(r,i) c 2 „ „ „, . „„. , 
—V x V x E(r, t) + 2E(r, t) 



uj p2 dt 



p2 



(5.11) 



_ 4 n^Ei0 E(r) t j _ 4^2_ E(r; t) x B s (rj t) = Q 



"0 



VI. NON-LINEAR COUPLING FIELDS EQUATION 



Now wc study the low- frequency field equation l|3.4[) . Differentiating with respect to section v in braces in 
and taking account to delta function in Ea. il3.4[) . yields 



S. 



2(t) 



1 gf n o ( ^ 



4m|u; 2 2 w \^2kv T 2 



(k-eLJ-efc, (k-<)]. 



Then Eq. I|3.4[) becomes 



(l 2 2 2 A \ -ciTS 

(fee — u> e k ) E fc 



iq 2 ( lu 



m 2 \ \[2kv T 2 



Jkx (e^ (+) x E£ _) )<f(fc - fci - fc 2 ) dfcidfc 2 . 



According to Ea. (|4.t)|) . 



we have 



(fc 2 c 2 - ^ 2 4) Bf 



e"^Bf = -k x E^ 5 , 



k X 



k X 



E^ (+) x E T k ^ 



^ (5 (fc — fci — fc 2 ) dkidk2 



For very low frequency fields, «ti 3> w/fc ^cC wt2, one has 



To 



fc 2 l> 



U) TT U) ITT U) 

- i\hr~, ~ —i\ 



2 kVTa 



2 kVTa 



pl 



it 2 '. 



Tl 



fc 2 4 2 



P 2 _. /7tw P 2 a; 



2 w kvTo 



I IT U. 
-l\ — 



2 UjkVT2 



and the following condition obviously is satisfied 



/ kvT2 i 
Up2 > I ) KC. 



Then Eq. I|(j.4[) become 



wB €e-* = -^k x 



m 2 o; 2 2 



k x / (e^ (+) x e£ 2 (_) )<5 (fc - fci - fc 2 ) dfcidfc 2 
Hence the coordinate representation of Eq. (|6.7|l is reduced to the following equation 



|-B s (r, t) = x V x [E(r, t) x E*(r, t)] , 

c/t m 2 uj p 2 

where we have chose = f for getting real magnetic fields with low-frequency. Through the substitutions 



E(£,r) 



3 fT2/w p2 ' 

4V3E(r,i) 



2 mi c 
T=-fj,u p2 t, (i= — , a = — 

3 m 2 -j 1 't2 

B(f,r) = r ^-B-(r,t), n=§£, 
A^m 2 cuj p 2 4 n 
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we can now 



write Eas. l|5.7|) . I|5.11|l and 16.8|l in the form 



n^,r) = m,r)\ 2 , (6.10) 



^ E ^' T) - «V x V x E(C, r) + |-E(£, r) - n(£, r)E(£, r) - iE(£, r) x B(£, r) = 0, (6.11) 

|-B(e, r)^)vxVx [E(e, r) x E*(£, r)] . (6.12) 
ar 6 

Therefore it may be seen that the gravitoelectric and gravitomagnetic fields, and self-generated gravito-magnetic 
fields with very low-frequency are completely determined by the closed Eas. (|6.10[l - H6.12l) . 

VII. NUMERICAL INTEGRAL 

It is well known that the gravitational system can emit energy in quadrupole radiation. For example, consider a 
neutron star with mass M = Mq, its radius, rotation inertia, rotational period and eccentricity are as follows: 

R~10km, M~M , J- 10 45 , P~ 0.033, e - 10~ 4 . 
As a result, due to the quadrupole radiation, the energy loss for the system is[l8| 

L w 10 45 e 2 ~ 10 37 (erg-s- 1 ). 
Furthermore, one can estimate the radiation flux at r = D = r{ = 10 x lOfcm = 10 (cm), 

F ~ L/D 2 ~ L/rl ~ 10 23 {erg ■ cm- 2 • s" 1 ), 
where is the inner radius of accretion desk around the compact star , at which the corresponding parameters are 

T a =T 1 =T 2 = 10 7 K, Pl = n Q m p = p 2 /W = A{g ■ cm,- 3 ). 
Then we get immediately the estimated values of the field in the gravitational wave 

|E(r,t)[ 2 |E(r,t)| 2 „„ , nl2 , _ 3 . 

r ~ c , ~ b jc ~ 3.3 x 10 [erg ■ cm ), 

on OTT 



i.e. 



E(r,t)| 2 _- r _ 3.3 xlO 14 _ 3.3 xlO 14 _ in _ 3 _ 



and 



hence 



= W = = = 10 

87rnoTo 1.38 x 10- 16 n T 1.38 x 2.4 x 10 16 



|E(r,t)|* =0 = — W- 

/I 



|E(e,r)|^ =0 = 3.84, a = = 3.8 x 10 6 , (TT r=0 = 10~ 3 ). 
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We have solved numerically Eas. H6.10|) . I|6.11|l and (|6.12() in two dimensions with three field components using FFT. 
The initial condition with a periodic boundary condition is given as |l6tll^| 



. l2 „ . 27TW X w , ^ Ho 27TW X , X , . 

E (?! r ) r=n = -Eosm secft(-— Xe^ + e z ) - E cos tanh — sec/i— e,, (7.1) 

2/o £ o 27rL 2/0 L L 



with 



L = 2 x 10 3 , y = 5 x 10 6 . 

The distribution of intial gravitoelectric field is shown in Fig.l, the dynamic evolution behavior for gravitoelec- 

tric fields and self-generated gravitomagnetic fields with the very low frequency is shown in Figs. 2-7 and Figs. 8-13. 
Quantites in Figs. 1-13 are dimensionless. The relations to dimensional ones are 

IT 

r = 3 x 10 6 J— x Fig = 1.5 x Vfx Fig , (7.2) 
V 

( = 1.6x 10 4 -J= = 2.5 x 10 3 t, (7.3) 



= 3.6 x 10- 20 (n T ) |E| 2 Fig = 8.6 x 10 11 |E| 2 Fig , (7.4) 



B = 5.6 x WB Fig . (7.5) 

Figs. 2-7 give the collapse development of gravitoelectric field and Figs. 8-13 the collapse development of self-generated 
gravitomagnetic field with very low frequency. While r — 1.65, the strengths of the self-generated gravitomagnetic 
fields are 



|E|L X - 2748.75 , IB]^ = 8.64 x 10" 
Using Eqs.(jZ3 ~ 17.5|) . while t — 3.3 x 10 3 s we have 

|E(r,i)|Lx = 5.94x W 16 erg/cm 3 

and 



B s = 5.1 x I0 5 (erg/cm 3 ) 

[ 

case the expansion Ea. (|2.7|) is no longer valid 



When t > 1.65, the collapses rapidly and leads to a very strong field, i.e. W = |E(r, t)\ 2 /87rnoTo > 1 and in this 



VIII. CONCLUSION REMARKS 

From the above study, we arrive at the following conclusions: If there is very intense gravitational radiation with 
high frequency near by a source, the interactions of the wave decay and fusion and wave-particle develop, and nonlinear 
(matter) currents with low frequency are induced by the interactions, leading to increasing of local matter density 
and excitation of a very low-frequency gravitomagnetic field. 

It is shown that the dynamic behavior and configuration for the GEM fields, the perturbed density and the gravit- 
omagnetic field with very low-frequency are determined by Eas- H6.10|) - H6.12fl 
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We investigate the numerical solution of Eas. H6.10|) - (|6.12() . It is shown that the fields, involving the GEM fields, 
the perturbed density field and self-generated gravitomagnetic field with very low frequency, may collapse, so that the 
gravitational waves are amplified by a factor of 10 3 . In other words, due to self-condensing, a stronger GME fields 
could be produced; and they could appear as the gravitational waves with high energy reaching on Earth. In this 
case, Weber's results, perhaps, are acceptable. 

Meanwhile, the increasing perturbed matter density and gravitomagnetic field with very low-frequency in the local 
region, where the gravitational waves get through, are in favor of the formation of a new rotating object.// 
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FIG. 2: (a) r = 0.05, |E|l ax = 3.895 



FIG. 3: (b) r = 0.5, |E|L X = 3.418 



FIG. 4: (c) r = 1.4, jE^ = 5.213 




FIG. 5: (d) r = 1.55, |E|^ x = 12.09 



FIG. 6: (e) r = 1.6, |E|^ = 30.55 




FIG. 7: (/) r = 1.65, |E|L X = 2749 



FIG. 8: (a) r = 0.05, |B|^ = 1.171 x 10" 




FIG. 9: (6) t = 0.5, |B|^ ax = 6.938 x 10~ 



FIG. 10: (c) r = 1.4, |B|^ = 5.705 x 10" 




FIG. 11: (d) r = 1.55, |B|L X = 1-154 x 10" 



FIG. 12: (e) r = 1.6, |B|i ax = 1.405 x 10" 




FIG. 13: (/) r = 1.65, |B|^ ax = 8.641 x 10 



